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Abstract 

In 1962 Dyson used a physically based, macroscopic argument to deduce the first two 
terms of the large spacing asymptotic expansion of the gap probability for the bulk 
state of random matrix ensembles with symmetry parameter f3. In the ensuing years, 
the question of asymptotic expansions of spacing distributions in random matrix 
theory has shown itself to have a rich mathematical content. As well as presenting 
the main known formulas, we give an account of the mathematical methods used 
for their proofs, and provide some new formulas. We also provide a high precision 
numerical computation of one of the spacing probabilities to illustrate the accuracy 
of the corresponding asymptotics. 

1 Introduction 

Random matrices were introduced in physics by Wigner in the 1950's jl9]. Wigner's 
original hypothesis was that the statistical properties of energy levels of complex nuclei 
could be reproduced by considering an ensemble of systems rather than a single system in 
which all interactions are completely described. This allowed for an entirely mathematical 
approach where statistical properties of the spectrum of an ensemble of random matrices 
were considered. But coming from physics, the aim was to use mathematics to compute 
experimentally measurable statistical quantities, and to compare against the data. 

One viewpoint on a real spectrum from a random matrix is as a point process on 
the real line. As such, perhaps the most natural statistical characterisation is that of 
the distribution of the eigenvalue spacing. This choice of statistic becomes even more 
compelling when one considers that in many cases of interest, eigenvalue spectra can be 
'unfolded'. This means that unlike many statistical mechanical systems, the density is not 
an independent control variable, but rather fixes the length scale only. Unfolding then 
is scaling the eigenvalues in the bulk of the spectrum so that the mean density is unity. 
It is indeed the bulk spacing distribution for the Gaussian orthogonal ensemble of real 
symmetric matrices — albeit in an approximate form known as the Wigner surmize (see 
e.g. jH]) — which was compared against the empirical spacing distribution for the energy 
level of highly excited nuclei (again see [44] , and references therein). 
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Fixing length scales at the edge of the spectrum practical exercise, a more 

difficult task. In addition to the bulk, we will have interest in the soft and hard spectrum 
edges when the eigenvalue spectrum exhibits a square root profile and inverse square root 
profile respectively. To specify realisations of the bulk and edge regions of the eigenvalue 
spectrum, we recall (see e.g. [33]) that the so-called classical random matrix ensembles 
have their eigenvalue probability density functions (PDFs) of the form 

1 N 

1=1 l<j<k<N 

with /3 corresponding to the underlying global symmetry (/3 = 1,2 or 4 for invariance 
under orthogonal, unitary or symplectic unitary transformations respectively); C denotes 
the normalization. This is extended to general > 0, giving the /3-ensembles [2T] as 
specified by the eigenvalue PDF fll.il) . to be denoted ME N ^(g(X)). In particular the choice 
g(\) = e~ /3A / 2 defines the Gaussian /3-ensemble and the choice g(X) = A /3a / 2 e~ /3A//2 , (A > 0), 
defines the Laguerre /3-ensemble. 

The bulk state can be realized by scaling A/ t— > xi/y/2N in the Gaussian /3-ensemble. 
The soft edge is realized by the scalings A; H- y/2N + ^^ 1/6 and A; (-)■ 4iV + 2\[2x\ in the 
Gaussian and Laguerre /3-ensembles respectively [29] . Only the Laguerre /3-ensemble has a 
hard edge, as it requires the eigenvalue density to be strictly zero on one side; it is realized 
by the scaling \i i— > f^. In all cases the limit N — > oo needs to be taken after the scaling. At 
an edge, the spacing between consecutive eigenvalues is not the natural observable. Instead, 
it is most natural to measure the distribution of the largest, second largest etc. eigenvalue 
(or smallest, second smallest etc.). It is well known, and easy to verify, that all these 
quantities can be expressed in terms of the (conditional) gap probabilities E^\n; J) for 
there being exactly n eigenvalues in the interval J, for the scaled state (■) = bulk, soft or 
hard indexed by 0. In the case of the hard edge, the probability depends on the exponent 
f3a/2 in the Laguerre weight A^ //2 e~ /3A / 2 , so we write i^ ard (n; J; a/3/2). 

Our interest in this review is on the asymptotic form of spacing distributions in the 
bulk, and of the distribution of large and small eigenvalues at the edge. This is a topic 
which (in the bulk case) occupied the attention of Dyson in one of the pioneering papers 
on random matrix theory in the early 1960's [22], and is still being written on as we stand 
today some 50 years later. We are seeking to catalogue both the results, and the methods 
which underlie them, and also to contribute some new formulas. Section [2] deals with re- 
sults founded on Dyson's heuristic physical hypothesis; these are in the form of conjectures. 
The various mathematical techniques which can both prove, and build on these asymp- 
totic expressions, are covered in Section 3. A numerical illustration of the accuracy of the 
asymptotic form is given in Section 4, as is a discussion of asymptotic results for the gap 
probability in the case that each eigenvalue is independently deleted with probability (!—£)• 
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2 Macroscopic heuristics 



2.1 Zero eigenvalues in the gap 

The eigenvalue PDF (II. ip can be interpreted as the Boltzmann factor of a classical log- 
gas system interacting at inverse temperature /3. The particles repel via the logarithmic 
potential and are subject to a one body potential with Boltzmann factor g(X) = e~^ v( - x \ 
This interpretation led Dyson [22] to hypothesize an ansatz for the asymptotic form of the 
gap probability Ep(Q; (— a, a); C/3En), where C(3En denotes Dyson's circular ensembles 
(see e.g. [33j Ch. 2]) of random unitary matrices (all eigenvalues are therefore on the unit 
circle; the interval (—a, a) refers to a sector of the circumference specified by its angles), 

Ep(0;(-a,a);C/3E N ) ~ e^ 5F . (2.2) 

N— >oo 

Here and below the symbol ~ is used to denote that the RHS gives leading terms, up to 
some order to be further specified, of the asymptotic expansion of the LHS. In ( 12. 2ft 5F 
is the energy cost of conditioning the equilibrium particle density so that P(i)(9) = for 
9 G (—a, a). This energy cost consists of an electrostatic energy 

Vi = ~\ f V r (p(i)(0i) - N/2n) (p (1) (0 2 ) - N/2n) log |e^ - e ie *\d9 1 d0 2 (2.3) 



and an entropy term 



1 P ^ 



V4 = (j-§;y i°s(f§) (2.4) 



The density is chosen to minimize V\ and then V\ and V2 evaluated, and we have 

5F=(V 1 + V 2 ). (2.5) 

Proposition 1. (Dyson 1962) With the requirements that p{\){9) = for 9 G (—a, a) and 
Jq 7 " p(i)(9)d9 = N, V\ is minimized by 

p a) (9) = v / 7 2.6 

( H; 2tt ^/sin 2 (e/2)-sin 2 (a/2) 1 ^ 



We then have 
PV 1 = 



Vlogcos(|), PV 2 = (l-0iVlog(sec(|)+tan(|)). (2.7) 



We remark that explicit calculations in [22] showed that requiring p{\){9) to minimize 
V1 + V2 (rather than Vi) results in a correction to PV 2 which for large N is of order \og(Na), 
indicating that the asymptotic expansion (12.21) will not correctly give terms of this order. 
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Substituting (12. 7p in (12. 5p . and substituting the result in (12. 2 p gives a large deviation 
formula, telling us (as a conjecture) the probability of there being no eigenvalues in the 
interval (—a, a). This probability decays as a Gaussian in N. An 0(1) expression should 
result from choosing the excluded interval as (—tts/N, tts/N), as then there are 0(1) eigen- 
values in the gap. Replacing a by na/2 in (12. 7p . then taking iV — >• oo (this is a double 
scaling limit) gives the prediction 

lim Ee(0; (-tts/N, tts/N); C(3E n ) ~ e -P{«>)'n*+W2-V*>/*, ( 2 .8) 

N— >oo 

Dyson was well aware that the ~ symbol should be interpreted as agreeing in the large s 
asymptotic expansion to the order given. But the LHS is the definition of E^ nlk (0; (— s/2, s/2)), 
thus providing the following conjecture. 

Conjecture 1. (Dyson 1962) We have 

< llk (0;(0,s)) ~ e -/3(-) 2 /i6+(/V2-iW2. (29) 

s— >oo 

As remarked above, Dyson [22] carried through the details of the minimization of 
V\ + V2, resulting in a logarithmic correction to the exponent of the RHS of (12. 9ft : ((1 — 
/5/2) 2 / (2/3)) logs. However, this was later put in doubt by des Cloizeaux and Mehta 
[4"5] who using a method based on eigenvalues (see Section 13.31 below) obtained —1/8, 
—1/4 and —1/8 for the prefactor of logs for /3 = 1, 2 and 4 respectively. In 1976 [23] . 
Dyson himself used inverse scattering methods applied to the Fredholm determinant form 
of E^ nlk (0; (0, s)) (see Section l3TTj) to also give the prediction —1/8 for the prefactor in the 
case (3 = 1. In fact the correct extension of (12.91) for general /?, as proved for the Gaussian 
j3 ensemble, is [60] 

£j ulk (0; (0, s)) ~ exp (-/3(tts 2 )/16 + (/3/2 - l)vrs/2 + (l/4)(/3/2 + 2//3 - 3) logs + 0(1)). 

(2.10) 

Its derivation will be reviewed in Section 13.21 

The ansatz (12.21) was applied to the gap probability at the hard edge of the Laguerre 
ensemble by Chen and Manning in 1994 [11]. They considered the probability of there 
being no eigenvalues in an interval (0,t). 

Proposition 2. (Chen and Manning 1994) F° r the Laguerre ensemble specified by U.l\) 
with g(X) = X a e~ x , with the eigenvalues constrained to the interval (t,b), with t > given, 
the minimizing solution for the level density pm (x) is 

*>M-55n/i5Hi/D- (2 ' u) 

Normalization of the density requires that b is related to N by 




(2.12) 
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Using (12.111) appropriate analogues of (12. 3p and (12. 4p were computed (see also [T2]). 
thus giving a prediction for the large N form of ^(O; (0,i); ME A r i/ g(A a e~ A )). This is expo- 
nentially small in N. But with t = s/ (4iV), the number of eigenvalues in (0, t) will be 0(1). 
With the resulting expression interpreted as the large s asymptotic form of .Eg ard (0; (0, s); a) 
(s must be scaled s i— > (/3/2) 2 s to account for the latter being defined as the large N limit 
of ^(0; (0, s/(4iV)); ME JV , /3 (A a e- /3A / 2 ))), the following conjecture was obtained. 



Conjecture 2. (Chen and Manning 1994) We have 

4 ard (0; (0, s); a) ^exp (-^ + a^7s - £ log s + (l - £)± log s) . (2.13) 



Historically, (12.13p had already been proved for a G Z> and 2//3 G Z >0 in [30] before 
the work of [11]. Moreover, the work [30], which was based on a-dimensional integral forms 
for i?l ard (0; (0, s); a), gave the explicit form of the constant term in the extension of ( 12 . 1 3[) 
to next order (see Section I3T41) . 

The first application of the log-gas ansatz (12.21) at the soft edge was due to Dean and 
Majumdar [131 US] . 

Proposition 3. (Dean and Majumdar 2006) Consider the Gaussian f3 -ensemble ME^^{e~^ Nx2 
Suppose the eigenvalues are confined to the interval (—b,t) where t < 1 and b > is deter- 
mined by charge neutrality. The corresponding density is given by 

, . 2N /l-t + x\V*/l 



tv V t-x J \2 

where I := b + 1 = §(t + Vt 2 + 3). 

Only the corresponding form of V\ was computed, and this gave the large deviation 
formula 

^(0;(t,oo);ME^(e-^ 2 )) 

'2t 2 t A 5 , = 1 , , = 1 t + 



N-too 



exp ( - (3N 2 { — - — - -tV3 + ¥ - -t 3 V3T¥ - - log 
P V ' \ 3 27 18 27 2 f 



(2.14) 



and correspondingly, upon the appropriate soft edge scaling y2N{l — t) = — s/(y/2N 1 / 6 ), 
the asymptotic formula 

Ef\0;(s,oo)) ~ e -P\s\ 3 /2\ (2.15) 

S— >— OO 



This latter prediction was already implied by earlier work [29], [5 
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2.2 Loop equations 



In 2011, Borot, Eynard, Majumdar and Nadal [8] gave an alternative heuristic formalism 
to the Dyson log-gas ansatz, for purposes of computing the soft edge gap probability. This 
in based on the so-called loop equations associated with the large N form of the multiple 
integral definition of the latter. The approach allows for the Dyson ansatz (12.21) to be 
extended to include higher order terms; in practice two new terms are computed — one is 
termed the Polyakov anomaly, and the following result is obtained. 

Conjecture 3. (Borot, Eynard, Majumdar and Nadal (2011)) We have 
£J° ft (0;( S ,oo)) s _~ ixj exp(-/3J|£+ 

m^A^, + m)Hm-s log |s| + logT „ ft + 0(|s| _ 3/2) ^ (2 16) 

where 

log , r = (E _ g«, fl/2) + log2 _ !2iM _ » + Km ( , 17) 

with Kp the constant term in the large N expansion of F(N + 1) := ^^L 1 logT(l + j/3/2). 
(Note that in J2J/ what we call (3/2 is written as (3.) 

In [8], for (3 rational, Kp was evaluated in terms of the Barnes G-function, while for 
general (3 > it was shown 

log(27r) , f3 ( 1 x 7 f 00 1 ( s s s\ 

K ^ = ^ + 2{-2- C{ - 1) ) + 6i3 + J ^3l(^T- 1 + 2-12J dS ' 

where 7 denotes Euler's constant. In fact K^j 2 can be expressed in terms of the so-called 
Stirling modular form p 2 (l,r), which from a computational viewpoint can be defined by 
the infinite product [51] 



00 Q(nr) 

p 2 (l,r) = (2 7 r) 3 /V^/4 + (. + iA)/i2 e PH "Q * 



, r(l + nr) 

n=l v ' 



where 



P(t) = -j^ - ^ + rC'(-l), Q(x) = (- + x) log 2 - x + log ^Phx 



12t 12 3 v " " v ' V2 / 12x 

The quantity ^(Ijt) is fundamental to the theory of the Barnes double gamma func- 
tion Vi{z\ l,r) [2], the latter being related to the usual gamma function through the two 
functional equations 

1 _r-/-i/a T (z/r) 1 _ 1 T(z) 



r 2 (* + l;l,T) r 2 (2;l,r)' r 2 (z + r;l,r) ^T 2 (z;l,r) 

and furthermore is normalized by requiring lini^o zT 2 (z; 1, r) = 1. 



Proposition 4. Let r — 2/(3 and specify F(N + 1) and Kp/ 2 as in Conjecture^ We have 

F(N + 1) = (2 7 r) jV / 2 r-( jv2 - jV ( 1 -))/ 2 - r(iV)r(1 + ^ /r) , (2.19) 

T 2 {N; l,r) 

«i/r = -^logr + log27r-logp 2 (l,r), (2.20) 



TOi/i i/ie latter equation substituted into (2.17) giving 

logr| oft = - |((/3/2) + (2//3)) log2 + _ i ogp2 (i, 2 /fj). (2.21) 

The equation (I2.19p has appeared in the recent work [ID] ; it follows immediately by 
characterizing ,F(JV + 1) as a first order recurrence, and using (I2.18p . The formula for 
K\i T then follows by extracting the term independent of N in the corresponding asymp- 
totic expansion. Here one uses the fact that for logr 2 (iV; l,r) this is logp 2 (l,r) [50] . A 
consequence of (I2.2ip is that 

log = _ log « = m±m^ log , ^ ^ i log ^ (2 . 22) 



r. 



2//3 



p 2 (l,/3/2) 8 b V2; 2 b 2 



where the final equality follows from the inversion formula for the Stirling modular form 
[4TJ1 Prop. 7(iv)]. Using this in ( 12161) gives that (cf. P Eq. (6.2)]) 

^ oft (0; ( S) oo)) ^ (|) V2 ^(0; ((f ) 2/3 ^, oo)), (2.23) 
where E% oit refers to the RHS of ( 12.161) with |s| 3 / 2 replaced by — |s| 3 / 2 . 



2.3 Conditioning n eigenvalues in the gap 

In 1995, Dyson [21], and independently Fogler and Shklovskii [27], further developed the 
log-gas argument by the consideration of the setting that the gap (—t,t) is required to 
contain exactly n eigenvalues, with C n <C i. Moreover, a change of viewpoint was 
introduced: the log-gas was taken to be infinite in extent, with the bulk state characterized 
by a uniform density, normalized to unity. The n eigenvalues are modelled as a continuous 
conductive fluid occupying the interval (—b,b) C (—t,t). The electrostatic potential in 
this region must therefore be equal to a constant —v say, v > 0, with the potential in the 
other conducting region M\(—t,t) taken to be zero. The explicit form of the density was 
determined, and this substituted in the appropriate modification of (12.31) and (12 .4p gave 
after some calculation the simple results 

Vi = -y + y(t 2 -& 2 ), V 2 = v. (2.24) 

The end point b is determined by n via a certain elliptic integral, and similarly v in terms 
of an elliptic integral of modulus b/t. Expansion of these quantities for t — > oo, and 
substitution in (12. 2ft provides a generalization of Conjecture 12. 10H . 
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Conjecture 4. (Dyson [24)1 , Fogler and Shklovskii [27\l (1995)) For 0<n<Cs we have 



log£j ulk (n; {0, 8 )) 



— »oo 



-(3- 



7TS 



16 



+ 



n / /3 /3n 
2 V ~ 2 ~ ~2~ 



if£ + --3 

4v2 + /3 d 



logs (2.25) 



(Tiere we nave added the n = contribution to the term logs as implied by A2.10\) — we 
t/ien expect A2.25\) to hold for < n <C s; i/ws zs noi a consequence of the calculations in 

Only very recently has this infinite log-gas formalism been applied to predict the asymp- 
totic forms of the conditioned gap probabilities at the hard and soft edges [21]. Since the 
system is (semi-) infinite, this relies on characterizing these edges in terms of the respective 
background densities: a/x/7t for the soft edge, and 1/(27Ta/x) for the hard edge. In both 
cases the coordinates are chosen so that the edge occurs at x = 0. It was found in [37] 
that applying the ansatz (12. 2 j) with SF given by ( 12. 5 j) in this setting to the n = case gave 
results inconsistent with both (I2.13P and its soft edge analogue in the second order term. 
Thus the ansatz (12. 2p with SF given by (12. 5p is incorrect in the infinite log-gas formalism 
applied to the hard and soft edges. On the other hand, it was observed that replacing Vi 
by the potential drop v in going from the region containing the infinite mobile log-charges, 
to the region containing the n charges — which according to ( I2.24p is an identity for the 
bulk — restores the correct value for these terms. Making this replacement for general n 
then gives the following predictions. 

Conjecture 5. (Forrester and Witte [S 7 ^) We have, for <^ n ^ \s\ (or more strongly 
< n< \s\), 



log£j ard (n;(0,s);/3a/2) ~ -pfi-y/jt 

1 s-s>oo lg \ 



a 

n H — 
2 



+ 



n 2 na a(a — 1) a 

\ h — H 

2 2 4 2/3 



logs 1 / 2 , (2.26) 



and 



\ogE^(n; (s,oo)) 



s— >— oo 



P\s\ 3 2^2, 



21 



P 2 
—n 

2 



1 3/2 



P 



H» + 5 HK 



log|sr 3/4 . (2.27) 



(As for \2.25\) . the results coming from the log-gas calculation have, in the case of the 
logarithmic term, been supplemented by knowledge of the asymptotic expansion at that 
order for n = 0.) 
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We remark that a check on (12.251) . (I2.26P and (I2.27P is that they obey certain asymp- 
totic functional equations, implied by exact functional equations for spacing distributions 
obtained in [52] • For example, at the hard edge one requires 

4 ard (n; (0, s/sp); Pa/2) -^jf {\(3{n + 1) - 1; (0, s/h/p) 5 a - 2 + 4//3) , 

where sp is an arbitrary length scale that satisfies S4/p(f3/2) 2 = sp. This is indeed a 
property of (I2.26p . 

Precise asymptotic statements can also be made concerning the asymptotic form of 
Ep\n; J), for |J| — > oo and n ~ {nj), where nj ((nj)) denotes the number (expected 
number) of particles in J for the unconstrained system. Thus macroscopic heuristics ap- 
plied to this linear statistic (see e.g. [MJ §14.5.1]) predict that (nj — (nj)) / 'y/Vainj has a 
Gaussian distribution with zero mean and unit variance, and so suggesting the following 
result. 



Conjecture 6. For n « (nj), 



p y ' |j|-Kx> (2vrVarn J ) 1 / 2 v ' 



Moreover, for (■) = bulk, soft and hard we have 



2(-sf/ 2 s 1 ' 2 
(n(o, s )) ~ s, (n {Sj0o) ) ~ — , (n (0 , s )) ~ (2.29) 

s— >oo s— >— oo s->oo 7T 



and 



2 1 1 

Varn ((M ) ~ -^-rlogs, Varn (Si0o) ~ -— log |s| 3/2 , Varn (0 , s ) ~ -— logs 1/2 . 



(2.30) 

The results (I2.29P are immediate consequences of the corresponding asymptotic density 
profiles (recall the second sentence below Conjecture H]), while f)2.30p can be derived heuris- 
tically from knowledge of the asymptotic form of the two-point correlation function (see 
[3"3"| paragraph below (14.87)]). In the case of (•) = bulk, (12. 28 p . with the corresponding 
vaues of (ri(o,s)) an d Varnj as implied by (I2.29P and ( I2.30p . was derived in the context of 
the infinite log-gas formalism by Dyson [24] and by Fogler and Shklovskii [27] . 



3 Rigorous methods 

3.1 Toeplitz/Hankel asymptotics 

It is a fundamental result in random matrix theory (see e.g. [33] §9.1]) that in the scaled 
limit (■) equal to bulk, hard or soft, and /3 = 2 the probability of there being no eigenvalues 
in an interval J, may be written in terms of a determinant of a Fredholm integral operator 

4 ) (0;J) = det(l-/4 ) ), 
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where K^) is the integral operator on the interval J with well known sine, Bessel and Airy 
kernels (see e.g. [33] for the precise definitions). This is related to the fact that for /3 = 2 
the gap probabilities can be written in terms of either Toeplitz or Hankel determinants. 
For example, the Toeplitz determinant of a function f(6), integrable over the unit circle, 
is defined as 



/ -, n2n \n-l 

D n (f) := det — / e«*-Wf{0)dB , (3.31) 

\ 2n Jo J jtk=0 



and one has the well known formula 

E 2(0 ; (-a, a); CUE„) = D N (U U={ \ ~ 

In particular lim^oo -D n (/ 2s / n ) = det (I— K^ nlk ), allowing for a strategy whereby the s — > oo 
behaviour can be extracted from the asymptotics of the Toeplitz determinant. On the other 
hand the Toeplitz determinant has a representation in terms of quantities associated with 
orthonormal polynomials (f>k(z) = Xk% + . • . with weight f(9) on the unit circle; explicitly 
D n {f) = rifclo Xfc 2 - Krasovsky jl3] used a Riemann-Hilbert formulation to compute the 
large n form of j- ln£) n (/ M ), uniformly in fi, providing both a proof and refinement of (12. 9p 
in the case (3 = 2. 

Theorem 3.1. (Krasovksy 2004, Ehrhardt 2006) We have 

log £ 2 bulk (0; (0, S )) = -M! - I log(7r S /2) + 1 log 2 + 3C'(-1) + O (^j (3.32) 

where ((z) is the Riemann zeta- function. 

We remark that up to the constant term this result, deduced by Dyson [23] using a 
scaling argument from known Toeplitz determinant asymptotics, was first rigorously proved 
by Deift, Its and Zhou [IE]; also the proof of Ehrhardt [25] is operator theoretic, and does 
not make use of orthogonal polynomials. 

Analogous strategies can be used to analyze the hard and soft edges for /3 = 2, giving 
the following results, proving and extending ( 12 . 1 3[) and ( 12 . 1 5[) respectively. 

Theorem 3.2. (Deift, Its and Krasovsky 2008, and Deift, Krasovsky and Vasilevska 2010) 
We have, for s — > — oo ]5 



\ogEf\0; ( a , oo)) = - J|L _ I l og |*| + 1 In 2 + C'(-l) + 0(| S |" 3 / 2 ) (3.33) 



and for s — > oo J7?)/ 



log E 2 hard (0; (0, s); a) = - S - + a^7s - j log s + log ( G ^tJ ) + 0(^ 1/2 ) (3.34) 



where G(x) denotes the Barnes G-function. 
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An alternative proof of (I3.33P has been given by Baik, Buckingham and DiFranco in 
2008 PP, using the Painleve form of i?f oft (0; (s, oo)) ^7\. This method carries over to the 
cases (3 = 1 and 4, and in [I] the expansion (12.161) with 

l 0g rf « = _ll]2g2 CM) „ tt = _371og2 CM) (3 35) 

6i 4g 2 , 6 4 4g 2 v ; 

was obtained. These confirm the values implied by (I2.2ip . 

With regards to (I3.34p . as noted above, for a £ Z> it was first proved by Forrester [30] . 
More recently a proof of (I3.34p valid for \a\ < 1 was given by Ehrhardt [26]. Furthermore, 
let the next order (constant) term in the exponent of (12.131) be included by adding log r^ rd . 
We read off from flBTMj) that r^ rd = G(l + a) / (2n) a/2 . For a £ Z> a multiple integral form 
for £'] iard [36], and an identity [35] relating E\ wd to i?2 ard and ^ ard for general a > -1 
tells us that 

_ 2 -a(q+i/2) g(3/2)G(2a + 2) _ Q(a+1)/4 -i r a , 2 

G(a + 3/2)G(a + 2) T(a-i)/2,i 

In the case of bulk scaling, include a constant term by adding log r^ ulk to the exponent 
of (I2.10p with s replaced by s/n (thus the bulk density is now l/n). It follows from (I3.32p 
that logT2° ulk = | log 2 + 3£'(— 1). And inter-relations between the bulk gap probability for 
= 1 and 4 with = 2 quantities give that [3], 

bulk = 2 5/12 e (3/2)C'(-l) ; r bulk = 2 -29/24 e (3/2)C'(-l)_ (3.37) 

We observe that (13.371) is consistent with a relation analogous to (I2.2ip . 
Conjecture 7. Let p 2 (l,r) denote the Stirling modular form. We have 

log^/f = (3 - \{0/2 + 2//?)) log 2 + 3(i log27r - log 0,(1, 2//3)), (3.38) 



and consequently 



. P2 (l,2//3) 



E^frfrs/*))^-) ^(0;(0,| S /7r)), (3.39) 



where Eg ulk refers to the RHS of $2.10\) with s replaced by —s in the second term. 



3.2 Stochastic differential equations 

The Gaussian and Laguerre /3-ensembles, defined as eigenvalue PDFs below ( 11. lj) . admit 
realizations as real symmetric tridiagonal matrices |21j . In the scaled N — > 00 limit, this in 
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turn leads to explicit characterization of gap probabilities in terms of stochastic differential 
equations. The first result of this type was done for the soft edge, by Ramirez, Rider and 
Valko [52] • With N fixed, it relies on expressing the number of eigenvalues greater than [i 
as the number of sign changes of the shooting vector for the tridiagonal matrix. Similarly 
at the hard edge (5Tj. In the bulk, the shooting eigenvector must be parametrized in terms 
of the corresponding Priifer phase |42l [59]. The following results are obtained. 



Proposition 5. Let b t denote standard Brownian motion. At the soft edge, define a diffu- 
sion by the Ito process JM 



2 

dp(t) = -=db t + (A + 1 - p 2 (t))dt, p(0) = oo; 
VP 



at the hard edge with parameter /3(a+ l)/2 — 1 by |5 

dp{t) = db t + (j(a + i) - ^V\e~ m coshp(t))dt, p(0) = oo; 
and in the bulk by I6U$ 

dp(t) = db t + (\tuahp(t) - f Ae _/3t/4 coshp(t)W p(0) 
V 2 8 / 



oo. 



Let J = (0, s/2vr) for (■) = bulk; J = (0, s) for (•) = /iard; J = (s, oo) /or (•) = soft. We 
have 

£^(0; J) = Pr(p(t) > -oo, Vt GR + U {oo}). (3.40) 

The utility of these characterizations for the purpose of asymptotics is that they allow, 
via the Cameron-Martin-Girsanov formula, f)3.40p to be rewritten as the expectation of a 
functional of a transformed stochastic process. In contrast to (I3.40p . this functional allows 
for a systematic, rigorous s — > oo asymptotic analysis resulting in a proof of (I2.10p - 
giving in the process the correct form of the general /3 > 0, logs term, for the first time 
- and a proof of (12. 131) for general /3 > and a > — 1. At the soft edge only the leading 
asymptotic form (12.151) has been proved using this approach [52] . 

For the large iV limit of the circular /3-ensemble, the Priifer phase has been used to prove 
the analogue of the Gaussian fluctuation formula (12.281) . namely Ep(n, (—a, a); C/3Eat) ~ 
(l/(27rVarri(_ QjQ ,))) exp(— (n— Na/n) 2 / (2Var n^ a>a \)) , where Var nr- a>a ) ~ (l/(7r 2 ,5)) log 



3.3 Predholm determinant/eigenvalue forms for Ejj\n, J)/E^\o, J) 

With (•) denoting bulk, soft or hard, let E^\j; £) be the generating function for {E^\n; J)}, 
so that 

oo 
n=0 
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Generalising the Fredholm determinant expressions for E^(0; J) from Sect ion 13. 1[ one has 
that for = 2 

oo 

E ( 2 \j;0 = det(l -£!$>) = -£A,), (3.42) 

1=0 

where 1 > Ao > Ai > A2 > • • ■ > are the eigenvalues of Kj. Consequently 

E^jn; J) _ Xji ■ ■ ■ Xj n , . x 

It has been known since the work of Gaudin [3H] that associated with K^ 1 ^ is a commuting 
differential operator. Furthermore, the work of Fuchs [38] uses this, together with a WKB 
asymptotic analysis, to deduce the s — > 00 asymptotic form of Aj (j fixed). It was noted by 
Tracy and Widom [56] that the latter implies the term with (ji, j2, ■ ■ ■ , j n ) — (0,1,..., n— 1) 
dominates as t — > 00. These authors carried out a similar analysis in the soft and hard 
edge cases [53, [58], so arriving at the following result (stated as Prop. 9.6.6 in [33]). 

Proposition 6. Let G(x) denote the Barnes G-function. For n fixed 



£ 2 bulk (n; (0,3)) 



G(n + l)n- n/2 2- n2 - n (ns)- n2/2 e n7rs 



£ 2 bulk (0; (0,s)) ^00 

Ef'jn; (Q,a)) G(n + 1) ^ /to/ ^ ' 

E| oft (0;(0,s)) s~oo 7r n/ 2 2( 5 " 2 +™)/ 2 1 7 ' P ' " 

^ ard (n; (0,a)) G(a + n + l)G/(n + 1) _ n(2n+2a+1) _ w 2 /2 _ an/2 2rt%/i . . 

£ 2 hard (0;(0,s)) ^00 G(a + 1) ' l J 

In [331 §9.6.2], as t -»■ 00, £ bulk (n; t) and £ bulk (n; i) are related to £ 2 ard (-; •) for partic- 
ular choices of the parameters. The asymptotics of the latter are known as noted in the 
above proposition, allowing us to extend the first result in (I3.44p to f3 = 1 and 4 ([33], eqns 
(9.100) and (9.102)). 

Proposition 7. For n fixed and /3 = 1 and 4 we have 

E^ ulk (n; (0, s )) e^ s l 2 ( fl 



where 



c P,n ,_, Rn2/i + (R/9 ._ Un/9 1 + - , (3.45) 



£ bulk (0; (0, S)) P ' H ( 7rS )^ 2 /4+(/3/2-l)n/2 



Cl 



G(n/2 + l/2)g(n/2_+l) 7r _ n/22 _ w(n+1)/4 



G(l/2) 

_ G(n + 3/2)g(w + l) ^_ no _ 2n(n+1) 
C4 ' n " G(37^ 
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According to the first asymptotic formula in ( 13.441) . (I3.45P is, for a specific C2 in , valid 
too for (3 = 2. Furthermore the functional form 13.451 for general (3 > coincides with the 
log-gas prediction (12.271) . and thus validates the latter for (3 = 1, 2 and 4, and furthermore 
extends it by the evaluation of cp >n . 

We would like to extend Proposition [7J to the soft and hard edge cases. For this, let 
V~ for (■) = soft, hard and J = (0, oo), (0, 1) respectively denote the integral operators on 

J, dependent on a parameter s, with kernels Ai(x + y + s) and ^Jaiy/sxy)- Write 

E { ± \t;J) = det(lT yfivf) 

and define 

Results contained in [31], and further refined in [6], tell us that for s — > oo 

Er ft (2fc;( S ,oo)) = Ef\k;(s,oo)) + ... 
Ef t (2k + l;(s,oo)) = ~£! oft (;( S ,oo)) + ... 

El oit (k;(s, oo)) = £ 1 soft (2£; + l;(2 2/3 s,oo)) + ... (3.46) 

Here terms not written on the RHS are exponentially smaller (in s) than the given term. 
To proceed further requires a property of the eigenvalues of V~ which although supported 
by numerical computations, to our knowledge is yet to be proven. 

Conjecture 8. Let Vj(j = 0,1,2,...) denote the eigenvalues of V- , ordered so that 
\ v o\ < \ v i I < \ v 2 I < • • • ■ Then v^j > while i4j + i < f or eac ^ J = 0; 1; 2, ... . 

It is well known, and easy to verify (see e.g. [331 §9.6.1]) that (v^f = \V , where {Aj } 

are the eigenvalues of Kj\ This fact, together with Conjecture [8] and the analogue of 
(13131) relating to£| oft (£;; (s, oo)), tells us that for s oo 

Ef\2k- (s,oo)) 
£f oft (0; (s,oo)) 
El oft (2k + l; (s,oo)) 
£| oft (0; (s/2 2 / 3 ,oo)) 

Knowledge of the explicit asymptotic form of A^ from [57], together with the asymptotic 
form of El ott (0; (s, oo))/E s A ° ft (0; (s/2 2 / 3 , oo)) implied by (12361) and (13T35D then allows us to 
extend the second result of (13.441) to (3 = 1 and 4. 



1 

(l_ A s)(l-A^)...(l-A s 2fc _ 2 ) 
1 

(1-A!)(1-A|)...(1-A^ 1 ) 



+ 



+ . . . 



(3.47) 
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Proposition 8. (under the assumption of Conjecture^) 
We have 

Ef'jn; (s,oo)) 



Ef ft (n; (s/2 2 / 3 ,oo)) 
£f oft (0; (s/2 2 / 3 ,oo)) »-+-oo 

At the hard edge, formulas structurally identical to (I3.46P hold [31], [6], with the impor- 
tant qualification that the additional label need to specify the hard edge gap probabilities 
is (a — l)/2 on the LHS of the first two equations, and a + 1 on the LHS of the third 
equation; on the RHS's it is a, a and a — 1 respectively, and in the third equation s is 
scaled by 4 instead of 2 2 / 3 . The analogue of (13.471) then allows the analogue of Proposition 
[8] to be deduced. 

Proposition 9. (under the assumption of Conjecture^) 
We have 

E^ d (n; (0,s);(a-l)/2) (n _ 1+o)/2 A G((w + p)/2)G((w + g + o)/2) 

£^ d (0;(0, S );(a-l)/2) *~oc 1 J 11 G(p/2)G((p + a)/2) 

xs -(n 2 +n(a-I))/4 e ii^ 

£ 4 hajd (n; (0, g /4);a+l) ^ ^ ^^V^tt) 1 / 2 ^ hard (2n + 1; (0, s); (a - l)/2) 
^ ard (0;(0,s/4);a+l) *~oc 6 T((a+l)/2) E^ rd (0; (0, s); (a - l)/2) 

As a final point in this subsection, we remark that the Gaussian fluctuation formula 
(1 2 . 2 8 f) can be proved for /3 = 2, using only the determinantal structure (13.421) and the fact 
that Var n j — > oo [TBI 155] . 

3.4 Hard edge — generalized hypergeometric functions 

In the case that a £ Z + and general /3 > 0, the hard edge gap probability -E^ ard (0; (0, s); a) 
permits evaluation in terms of a generalized hypergeometric function based on Jack poly- 
nomials pjt\zi, . . . , zjy). The latter are labelled by a partition K4 > > • • • > ftjv > 
of non-negative integers, and depend on the parameter a. For a = 1 they are the Schur 
polynomials, while for a = 2 they are the zonal polynomials of mathematical statistics; 
their precise definition can be found in e.g. [33], §12.6]. Defining Ck(zi, . . . , zn) as pro- 
portional to P« (zi, . . . , zjv) with a specific proportionality depending on a and k [331 
eq. (13.1)], and the generalised Pochhammer symbol [33l eq. (12.46)], the generalized 
hypergeometric function p Fq** is specified by (see e.g. [33j §13.1]) 

p FW( fll> . . • , a p ; fc, . . . , b q - Xl , . . . , x m ) := £ o j^r^ ^'^. ■ ■ ■ > *"0 - ( 3 - 49 ) 



G(n/2 + l/2)G(n/2 + 1) 



xexp 



W 2 G(l/2) 
T V 2 



2-(5/8)n 2 +(l/8)n^_ s ^-(3/8)n 2 + (3/8)n 



^ e _^ ( _ s) 3/ 2 F 1 sott (2^ + l;( g ,oo)) 
£? oft (0;(s,oo)) 



(3. 
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Like their classical counterpart, these exhibit the confluence property 
lim p F (a) (ai, . . . , a p ; b x , . . . , b q ; — , . . . , — ) = p - 1 F^ a \a 1 , . . . , a p _i; h, ... , b q ; x x , . . . , x m ). 

a p ->oo * dp dp 

Using this in the case p = q = 1, together with an integral expression for 1F1 (331 §13.2.5] 
we can readily express the conditional gap probability £?& ard (n; (0, s); a) for a, (3 E Z> in 

5 /2 

terms of the generalized hypergeometric function qF{ , extending the n = result of [30j . 



Proposition 10. Let f3a/2,/3 E Z> . We have 

Ef rd (n;(0,s);(3a/2) = A p {n, a ) s ™+(/V2Mn+a-i) e -/? s /8 

d Vl ... dy n \\{l- yj f a l 2 J] l^-y/ 

^° j=l l<j<fc<n 

x Ff /2) (.; a + 2n; (s/4)^ 2 , • • • , (sy n /4f),(3.50) 

where 

Mn , a) = ^ ( P_y(PY^ w±mr (3 . 51) 

and in the argument of oFf the notation (u) r means u repeated r times. Furthermore, in 
the case n = 0, this remains valid for general (3 > 0. 

An integral representation of qF^^ allows for the rigorous determination of the x — > oo 
asymptotic expansion of qF{ \j, c + 2(m — 1)/ (3; (x) m ) for c,2//3 E Z + [30], implying the 
corresponding asymptotic form of £^ ard (0; (0, s)). 

Proposition 11. (Forrester 1994) Let 2/(3 E Z + , and a/3/2 = m E Z> . For s — > oo we 
have 



, x m(m+l)/2/9-m/4 , , , 

4 ard (0; (0, S ); m) = W ( ~ ) ^/sW^ j x + Q ( J_ ) | ? (3 52) 



w/iere 

m/2 



hard _ 9 (2//3-l)m. 



/ i \m/-z m 

d W^J/P). (3-53) 



We see that (13.521) is in agreement with the log-gas prediction (I2.13P for general a > 
—1,(3 > 0, and furthermore gives the explicit form of the constant in the asymptotic 
expansion (to use (13.531) for m ^ Z> and check for example (13.361) requires an appropriate 
rewrite of the product using (12.181) . 
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8 /2 

In the case 

/3 = 4of iT , an integral representation not available for general (3 shows 
that for s — )• oo and y±, ■ ■ ■ ,y n ~ 1, [S] 

oFf /2) (j c; ( S /4)/V 2 , ( syi /4)^, . . . , (W^) 

= oFi m (_; c; (s / ' ^^y^U^)/* ^ + ^J_^ j . (3 . 54) 

This, substituted in (I3.50P implies, as a conjecture, the extension of the asymptotic formula 
(12.271) to include the constant term. 

Conjecture 9. For fin E Z> 0; let 

hard , _ 2 -("+»)^ /gx n (a+ n-lW2 £(o + 2j/g) IIj=0 £(1 + (j + 1)1/2) 

n! 11 (27r)1/2 n^ra+a+a)^)- 

For s — ?• oo we have 
E h ^ d (n; (0,s);/3a/2) 



(n) exp ( - /3 { - + + ^) log ^ }) (l + O 



£^ d (0;(0,s);/3a/2) 'W^^'^^ ^\ v tJ,t 1 V 2 ' 2 7 6 // V Ks 1 / 2 

(3.55) 

We can check that (I3.55f) is consistent with the results of Proposition [9j 
3.5 Approach to unity of £^(0; J) for | J\ 

Generally the gap probability is given in terms of the /c-point correlation functions {p^} 
according to 

4° (0; J) — 1 — J ( x ) rfx + 2l / / p (2) ( x ' ^ rfxc ^ 

Thus the leading |J| — >■ asymptotic form of (0; J) is determined by the asymptotic 
form of p^ rd (x) for x 0, p s ff{x) for a; — > oo and p^ lk (x, y) for x, y — > (for (•) = bulk, 
P{i){x) = 1 and so gives no distinguishing information). The calculation of the first and 
third is elementary [28], [30], while direct calculation of p s ff(x) is only known for = 1, 
and /3 even [19J. Collecting these together, we have the following result. 



Proposition 12. Let 

A - j-fa+iya/gyia+i r(i + £/2) _ (r(/3/2 + i)) 3 

^ - 4 UV^ r(1 + a)r(1 + a + /3/2) > ^ - r(/3 + i)r(3/3/2 + 1) • 

For t 0, 

Fj ard (0; (0, t); a) = 1 - / s a ds + 0(t a+2 ) 

Jo 

Fj ulk (0; (0, *)) = 1 - i + \b p j [ ( Sl - s 2 f d Sl ds 2 + 0(^ +3 ), (3.56) 

^ jo Jo 
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while for t — > oo 



t(i + b/2) r°° e~ 2l3x:i/2 / 3 / r°° p- 2 ^ 3/2 /3 N 

<(0;(t,oo)) = l- l (4 ^; 2 / xm^ dX + (j t -X^T dX \ (3 - 5?) 

Two distinct derivations of (I3.57P for general /3 > are known, both involving use 
of a non-rigorous double scaling limit [3H |9]. In [20], the stochastic differential equation 
characterization (recall Section 13. 21) is used to give a rigorous proof for general f3 > 0, but 
without determining the prefactor of the integral. 



4 Other aspects 
4.1 Numerical results 

Bornemann [7J E] has given a detailed study of the numerical analysis relating to the 
precise numerical evaluation of spacing distributions for /3 = 1, 2 and 4, working from 
the Fredholm determinant forms. As an end product he has provided a suite of Matlab 
programs implementing the theoretical procedures. The implementation in Matlab, with 
the arithmetic done in the hardware, means that the tails of the spacing distributions 
cannot be computed: their numerical values written as decimals are typically smaller than 
10~ 15 , and so double precision arithmetic typically truncates significant nonzero digits, 
leading to unreliable results. But with there being numerous exact and conjectured results 
relating to the asymptotics of spacing distributions, there is much interest in implementing 
the theory of [TJ E] using an arbitrary precision package. As a start, we have done this for 
the Fredholm determinant form for i?2 Ulk (0; (0, s)) (in fact we have modified the procedure 
of [7J E] by using instead of Gauss- Legendre or Clenshaw-Curtis quadrature rules, the 
tanh-sinh quadrature rule (see e.g. [61])). As a result we are able to tabulate 

ig-(0;(0,,)) 

r[s) £r k (0;(0,s))' (4 - 58j 

where £ 2 b ' as (0; (0, s)) is the asymptotic form of ^ ulk (0; (0, s)) as given by ( I3.32p , extended to 
the next two terms: l/(8(irs) 2 ) + 5/(8(7rs) 4 ) (these follow from the Painleve transcendent 
characterization of i?2 Ulk (0; (0, s)) (see e.g. [331 §9-6.7])). The values in Table [1] clearly 
illustrate the accuracy of the asymptotic expansion, even for relatively small values of s. 



4.2 Diluted spectra 

For a general one-dimensional point process, the generating function (I3.4ip can also be 
interpreted as the probability that there are no eigenvalues in the interval J, given that 
each eigenvalue has independently been deleted with probability (1 — £). In this setting 
the | J | — > oo asymptotics can readily be deduced, by making use of a heuristic analysis 
based on (E2BD 0]. 
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s 




1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



1.0046735914726577 
0.9998383226940526 
0.9999753765440204 
0.9999961026171116 
0.9999991096965057 
0.9999997235559452 
0.9999998946139279 
0.9999999537746553 
0.9999999775313906 
0.9999999881794448 



Table 1: Tabulation of the ratio of the asymptotic to exact bulk gap probability for (3 = 2. 



Conjecture 10. For < £ < 1 we have 



where (rij) is given by Ii2.29\) for (■) = bulk, soft and hard. 

We see from fl2T29l and (E2D, ( 127131) . f l2~15D that as a function of s the decay exhibited 
by (I4.59P is proportional to the square root of the leading decay of Eg . A method to 
prove (14.591) for (3 = 2, making use of (I3.42p . has been given by Pastur and Shcherbina 
|48j . Alternatively, for this (3, (I4.59P can be verified by using known asymptotics of the 
Painleve transcendent evaluations, as done for (•) = soft in [5]. 

An interesting feature of the asymptotic expansion of the relevant Painleve transcen- 
dents with < £ < 1 is that they contain oscillatory terms, in contrast to their asymptotic 
expansion with £ = 1. It is indeed the case that oscillations can clearly be seen in plots of 
^Ef ft ((s, oo); £) with < £ < 1 |5j. Dyson [24] has combined Coulomb gas and Painleve 
theory to deduce the asymptotic form ^"^((O, s); £) when £ — > 1 and simultaneously 
s — > oo, which is shown to involve an elliptic theta function; for fixed £ the asymptotic 
expansion of the relevant Painleve transcendent [46] involves only trigonometric functions. 
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